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THE PATH RESISTANCE METHOD FOR BOUNDING THE SMALLEST NONTRIVIAL 

EIGENVALUE OF A LAPLACIAN* 

STEPHEN GUATTERY*, TOM LEIGHTON*, AND GARY L. MILLER 5 


Abstract. We introduce the path resistance method for lower bounds on the smallest nontrivial eigen- 
value of the Laplacian matrix of a graph. The method is based on viewing the graph in terms of electrical 
circuits; it uses clique embeddings to produce lower bounds on A 2 and star embeddings to produce lower 
bounds on the smallest Rayleigh quotient when there is a zero Dirichlet boundary condition. The method 
assigns priorities to the paths in the embedding; we show that, for an unweighted tree T, using uniform 
priorities for a clique embedding produces a lower bound on A 2 that is off by at most an 0(log diamcter(T)) 
factor. We show that the best bounds this method can produce for clique embeddings are the same as for a 
related method that uses clique embeddings and edge lengths to produce bounds. 

Key words. Laplacians, Laplacian spectra, graph eigenvalues and eigenvectors, graph embeddings 

Subject classification. Computer Science, Applied and Numerical Mathematics 

1. Introduction. In this paper we consider methods based on graph embeddings for estimating the 
smallest nontrivial eigenvalue of the Laplacian matrix representation of a graph. The Laplacian is one of 
many ways to view a graph as a matrix; it is defined as follows: Let G = ( V . , E) be an undirected graph with 
vertices v lf ...,v n . Then the Laplacian of G is an n x n matrix L such that 




! 


degr cc(vi) 

-1 

0 


if i = 3 
if (i,j) € E 
otherwise 


It is not hard to see that L is positive semidefinite (all eigenvalues are > 0). Since the row sums of L are all 
zero, the smallest eigenvalue is zero. If G is connected the second smallest eigenvalue A 2 is positive. 

The study of the connection between Laplacian spectra (particularly with respect to A 2 ) and properties 
of the associated graphs dates back to Fiedler’s work in the 1970’s (see, e.g., [Fie73] and [Fie75]). These 
properties have been used in graph algorithms, particularly algorithms for finding small separators [PSL90]. 

The Laplacian also has an important role in representing physical problems. It often occurs in finite 
difference, finite element, and control volume representations of problems involving elliptic partial differential 
equations. These problems often include a Dirichlet boundary condition that specifies that the values at a 
set of vertices are zero. To represent this condition in the Laplacian, the rows and columns corresponding to 
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the boundary vertices are deleted from the matrix. The resulting matrix is positive definite, and its smallest 
eigenvalue is the one of interest. 

Bounds on the smallest nonzero eigenvalues of both forms of Laplacian have other important applications. 
Since the matrices are symmetric, their extreme eigenvalues can be used in computing their conditionnum- 
bers which are used in the study of iterative linear system solvers to estimate rates of convergence [HY81J, 
and to analyze the quality of preeonditioners [Axe92, GMZ95], Bounds on A 2 are useful in the analysis 
of spectral partitioning, both because A 2 occurs in bounds on cut quality [Moh89], and because they can 
be used in isolating structural properties of the eigenvectors used in making the cuts [GM95, ST96], The 
eigenvalue A 2 has been related to expansion properties of graphs, and can be used in determining if a grap 

is an expander [AM85a, AI 086 ]. 

Related work involving graph embeddings has been used to bound the mixing time for random walks 
[JS89, SJ89, DS91, Sin92, Kah96]. In this case, the bound is on the second largest eigenvalue of the Markov 
chain'transition matrix. Mixing time bounds have been used in approximation algorithms for a number of 
problems (see, e.g., [DFK91]). 

In many of the preceding applications, it is necessary to show bounds for classes of graphs in order to 
state results in terms of asymptotic algorithm behavior. For example, in the analysis of spectral partitioning, 
[GM95] uses eigenvalue bounds on a family of bounded-degree graphs to prove facts about the structure of 
eigenvectors used in partitioning; [ST96] gives an upper bound on planar graph eigenvalues that can be 
applied in bounds on the cut quotient of the resulting cut. The embedding techniques we present arc 
well-suited to producing such general results, and can be used with known results about embeddings. We 
have used them to generate lower bounds for families of graphs where the critical path resistance can be 

parameterized in terms of the size of the family member. 

Note that the applications listed above involve both upper and lower bounds. Upper bounds on A 2 are 
usually given by exhibiting small Rayleigh quotients. If x ^ 0 is an n-vector its Rayleigh quotient is ■ In 
fact A 2 = min v _ 0 ■ Thus a Rayleigh quotient for a vector x such that E x = 0 S lves an u PP er bound 

on A 2 . As a b^utiful example, Spielman and Teng have proved that A 2 = 0(l/n) for bounded degree planar 
graphs by showing how to construct a vector with a small Rayleigh quotient for each such graph [ST96] . 

In the zero Dirichlet boundary case we are interested in the smallest eigenvalue A. Consider the Laplacian 
L with the rows and columns of the boundary vertices deleted. Then A = min x x t x • bl r,b ' s case tbe 
Rayleigh quotient has a continuous analog: Let Q be a compact domain in !R d with boundary T, e.g. the 
unit square in the plane, and let u ± 0 vary over all differentiable functions defined on Q such that u - 0 on 
T. Then the Rayleigh quotient is 


where Vu is the gradient of u. This lower bound on c is known as a Poincare inequality. This analog provides 
useful insight on the problem and helps with the proofs to follow. 

We start by describing a very simple yet powerful method for bounding A 2 from below (Figure 1). We 
show that the Path Resistance Method correctly bounds A 2 from below. Note that the quality of the lower 
bound improves with the quality of the path embedding since better embeddings decrease congestion. We 
also show that it is sometimes possible to improve the lower bound by assigning a priority to each path and 
allocating the edge conductances based on the priorities. The path resistance method as stated in the figure 
does not mention priorities; we refer to that version as the unit-priority case. 
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Procedure: The Path Resistance Method 


Method: 

1* Construct a path embedding of K n into G , i.e., 
construct a path between every pair of vertices of 
G. 

2. For each edge in G compute the congestion 

Qj = \{P\ eij e P}\ 

3. For each path P and each edge e tJ on P allocate 
a resistor of size to P. 

4. For each path P compute its resistance, i.e., 
J2 e%J eP c ij- Let r be the maximum resistance over 
all paths. 

5. return “ n/r < \ 2 

Fig. 1.1. The Uniform Priority Path Resistance Method 


The clique embedding version of the Path Resistance Method is closely related to another clique embed- 
ding technique presented in [Kah96] that assigns edge lengths and computes the sum of the lengths of all 
paths incident to each edge. We show that these methods are duals of one another, and that the best lower 
bounds computed by these methods are the same. 

We also present a version of the Path Resistance Method that is applicable in the zero Dirichlct boundary 
case. Instead of embedding a clique, we embed a star into the graph, with the boundary vertices consolidated 
into the vertex at the center of the star. 

Finally, we show that the clique embedding version of the Path Resistance Method with uniform priorities 
is robust when applied to trees. We produce upper and lower bounds to show that for a tree T using uni form 
priorities, our estimate is off by at most a factor of 0(log diameter(T)). It is somewhat surprising that such 
a simple way of setting the priorities gives such a good bound. It is open as to how well the uniform priority 
method does for general graphs. 

2. Previous Work. Prior work related to the applications of the lower bounds has been cited in the 
Introduction. The use of clique embeddings to bound eigenvalues arose in the analysis of mixing times for 
Markov chains by Jerrum and Sinclair [JS89] [SJ89], Further work in this direction was done by Diaco nis 
and Strook [DS91] and by Sinclair [Sin92], Kahale [Kah96] generalized this work in terms of methods 
assigning lengths to the graph edges, and showed that the best bound over all edge length assignments is 
the largest eigenvalue of the matrix T T r , where T is a matrix representing the path embedding ([Kah96] 
also cites unpublished work by Fill and Sokal in these directions). He also gave a semidefinite programming 
formulation for a model allowing fractional paths, and showed that the bound is off by at most a factor of 

log n. He proved this gap is tight; he also noted that the results can be applied to Laplacians with suitable 
modifications. 

3. Terminology, Notation, and Background Results. Wc assume that the reader is familiar with 
the basic definitions of graph theory (in particular, for undirected graphs), and with the basic definitions 
and results of matrix theory. A graph consists of a set of vertices V and a set of edges E\ we denote the 
vertices (respectively edges) of a particular graph G as V(G) (respectively E(G)) if there is any ambiguity 
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about which graph is referred to. When it is clear which graph we are referring to, we use n to denote \V\. 

We use the term path graph for a tree that has exactly two vertices of degree one. That is, a path 
graph is a graph consisting of exactly its maximal path. A star is a tree with exactly one vertex that is not 
a leaf. We call the non-leaf vertex the center of the star. 

3.1. Matrices and Matrix Notation. We use capital letters to represent matrices and bold lower- 
case letters for vectors. For a matrix A, a {j or [A] tJ represents the element in row i and column j; for the 
vector x, x t or [x], represents the i th entry in the vector. The notation x = 0 indicates that all entries of the 
vector x are zero; 1 indicates the vector that has 1 for every entry. For Laplacians, we index the eigenvalues 
of an n x n matrix in non-decreasing order: Ai represents the smallest eigenvalue, and A n the largest. We 
use the notation A 4 (L) (respectively A,(G')) to indicate the i th eigenvalue of matrix L (respectively of the 
Laplacian of graph G) if there is any ambiguity about which matrix (respectively graph) the eigenvalue 
belongs to. u, represents the eigenvector corresponding to A,. 

3.2. The Laplacian Matrix Representation of a Graph. A common matrix representation of 
graphs is the Laplacian. Let D be the matrix with da = degree(u t ) for v t e V{G), and all off-diagonal 
entries equal to zero. Let A be the adjacency matrix for G K = 1 iff («i.«i) € E(G), 0 otherwise). Then 
the Laplacian representation of G is the matrix L = D — A. 

The following are some useful facts about the Laplacian matrix: 

• The Laplacian is symmetric positive semidefinitc, so all its eigenvalues are greater than or equal to 

0 (see e.g. [AM85b]). 

• A graph G is connected if and only if 0 is a simple eigenvalue of its Laplacian (see e.g. [AM85b]). 

• The following characterization of A 2 holds (see e.g. [Fie73]): 

x T Lx 

(3-1) A2 = S^t- 

• For any vector x and Laplacian L of the graph G, we have (see e.g. [Fie75]): 

(3.2) x t Lx= ^2 ( Xi-Xj ) 2 

e(G) 

An edge- weighted graph has a real, nonzero weight ui x] associated with each edge («», vj) (we consider 
a zero edge weight to indicate the lack of an edge). Fiedler extended the notion of the Laplacian to graphs 
with positive edge weights [Fie75] ; he referred to this representation as the generalized Laplacian. Let 
w i3 be the (positive) weight of edge (i,j) in graph G. Then the entries of the generalized Laplacian L of 
G are defined as follows: l u is the sum of the weights of the edges incident to vertex v t ; for 1 ± j and 
(vi,Vj) G E(G), kj = -Wij, and l i3 = 0 otherwise. 

With the exception of equation (3.2), the properties listed above also apply to generalized Laplacians. 
A slightly modified version of (3.2) holds for the generalized Laplacian L: 

X^/.X — ^ ^ Wij {X{ Xj ) 

(v r ,v 3 )£E(G) 

We usually use L to denote the Laplacian of a graph G, and K to represent the Laplacian of the clique 
K (whether K refers to the graph or the Laplacian will be clear from context). 

3.3. Graph Embeddings. We start with the notion of embedding one graph into another. Let G and 
H be connected graphs such that the vertex set of 11 is a subset of the vertex set of G. An embedding 
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of H into G is a collection T of path subgraphs of G such that for each edge (v,,v y ) e E( //), F contains a 
simple path 7^ from v % to Vj in G. For full generality, we allow fractional paths in our embeddings: i.c., 
an edge (vj,Vj) € E(H) can be associated with a finite collection of simple paths from 17 to v y in G; each 
such path has a positive fractional weight associated with it such that the weights add up to 1. If a path 7 
includes edge e, we say that 7 is incident to e. 

In the unit-priority case, the congestion cong(e) of an edge e € E{G ) is the number of paths in the 
embedding incident to e (fractional paths using e contribute their fractional weight to this count). The 
resistance r(e) of an edge is its congestion divided by its weight: r(e) = 2 Mfel. The resistance of a path 
r( 7 ) is computed from the edge resistances in the same way that the resistance of a series of resistors would 
be computed in an electrical network. That is, r( 7 ) = £ e e 7 r(e). The largest path resistance r max is defined 
as r max = max 76 r rfr). The corresponding path is called the critical path. Note that r max is always 
defined with respect to an embedding; the particular embedding will be clear from the context, and thus we 
will not specify it in the notation. 

The definitions change slightly when we use path priorities. Let p tJ be the path priority for 7 ijt and let p 
be the vector of priorities. The congestion cong(e) of edge e is the sum of the priorities of e’s incident paths. 
Each edge has a resistance factor r(e), defined as the congestion of the edge divided by its weight. The 
resistance of edge e with respect to is The resistance r(7) of a path 7 is the sum of the resistance 
factors of the edges in the path divided by the path priority (i.c., r( 7lJ ) = -L r(e)). The maximum 

path resistance taken over all paths in F is denoted r max . The corresponding path is called the critical 
path. 

We note that using these modified definitions, the general method for computing lower bounds using 
priorities is clear: Compute r max according to the revised definitions. Then ; n < \ 2 . 

4 . Lower Bounds on A 2 for Generalized Laplacians. 

4.1. A Bound Based on a Clique Embedding. We now consider the case in which G is a graph for 
which we would like to bound A 2 of G’s (generalized) Laplacian from below. We will use a clique embedding 
to decompose the graph. The decomposed graph is then viewed in terms of an electrical circuit analogy. 

We start with an easy lemma that helps convey the nature of the bounding technique. Let L be the 
Laplacian of (edge-weighted) graph G, and let K be the Laplacian of the complete graph on n = \V(G)\ 
vertices. 

Lemma 4.1. Let a be a positive real number. If aL - K is positive semidefinite, then \ 2 {L) > 2 . 

Proof Assume that aL-K is positive semidefinite. Consider u[(ai- K)u 2 , where u 2 is the eigenvector 
corresponding to A 2 (L). Without loss of generality, assume that u 2 has been normalized to length 1 . Note 
that K = nl - J, where J is the n x n matrix with all entries equal to 1. Thus it is clear that every vector 
orthogonal to 1 is an eigenvector of K with eigenvalue n. Therefore 

u 2 (aL - K)u 2 = au 2 Lu 2 - u 2 Ku 2 = a\ 2 - n > 0 , 

where the last inequality holds by the assumption that aK — L is positive semidefinite. The lemma follows 

□ 

Thus, if we can demonstrate an a sufficiently large to make aL — K positive semidefinite, we can bound 
A 2 from below. To do this, we form an embedding (F, p) of the complete graph into G, then break aL-K into 
pieces consisting of edges in K and their corresponding paths in T. We use the following additional notation: 
E tJ denotes the Laplacian of the graph on V(G) that has only the single edge (i,j). L(j iy ) represents the 
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generalized Laplacian on V(G) whose edges are the edges of 7«j> and whose edge weights are set as follows: 
each path 7^ using an edge (g, h) gets a share of w gh proportional to cot ^( g ,h) • Thls works out to 

Wgh. Pij Pij 

cong {g, h) r gh 

It is easy to see that the weights assigned to all incident paths for ( 9 , h) sum to w gh . (If fractional paths are 
involved, then the matrix L{'Uj) is multiplied by the fractional coefficient of path 7 y •) 

We now state the theorem that gives us the technique for bounding A 2 : 

Theorem 4.2. For any generalized Laplacian L of a connected graph G with positive edge weights, and 
any clique embedding into G f 


T max 


Proof. 

Any embedding of the complete graph into G defines a natural decomposition of the Laplacian L in 
terms of the i 0 ’s: define the matrix L nonpat h as the Laplacian of the graph on V(G) that includes any edges 
of G not used in any path 7. It is easy to sec that 


L — ^7 L(^f) d* Lnonpath • 
7er 


Since L nonpat h is positive semidefinite, the following inequality holds for any positive a and real vector x: 
(4.1) 


x T ( £ (aLfa) - Ea) ]x<x r (aL- K) 

i<j 


X. 


Applying properties of linearity and rewriting gives 


(4.2) Y, ( axTL (^) x - xT £d x ) < xT ( aL - K ) x 

l<j 

Now consider the terms ox T L( 7tJ )x - x T E v x. By (3.2), x r E tJ x = (x t - x 3 ) 2 . Likewise, the term for the 
path can be written in terms of its edges and edge weights: 

ax T L( 7ij)x = a ^ — (x g ~ x h ) 2 

9h 


Since all the priorities and resistance factors are positive, the sum in the expression above can be interpreted 
as the energy dissipation of a series of resistances (the reciprocals of the ^-’s) given a set of potentials (the 
Xi’s) at the nodes between the resistances. It is well known (see e.g. [DS84]) that this quantity is minimized 
when the potentials at the internal nodes are consistent with Kirchoff’s law (this is easily seen through an 
application of the Cauchy-Schwarz inequality). At that minimum the sum is the square of the potential 
differences of the endpoints divided by the path resistance: 


y; —(x g -xh ) 2 > 


(Xj - Xj ) 2 

rilij) 


Thus, if a = r mttx , we have for every path 7 ij and every x that 


ax T L( 7y )x - x T E ij x > r m J X ‘ . - (x t - x 3 f > 0. 
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For this value of a, the left-hand sides of inequalities (4.1) and (4.2) are nonnegative. Thus r max L - K is 
positive semidefinitc, and the theorem holds by Lemma 4.1. 

If fractional paths are involved, then each L( 7^) is multiplied by the fractional coefficient, as is the 

corresponding E %} . The fractional coefficients cancel, so it is easy to sec that the theorem also holds in this 
case. □ 

We now compare the path resistance method with the edge length method described in [Kah96]. The 
method as applied to Laplacians is described as follows: 

• Specify a clique embedding for the graph, and assign each edge a positive length. 

• Compute the length of each path with respect to the edge lengths. 

• For each edge, compute the sum of the lengths of all incident paths divided by the length of that 
edge. Let p max be the maximum such value taken over all the edges. Then — — is a lower bound 

r \ Pmax 

for A2. 

The proof that this produces a lower bound can be done in a fashion similar to the proof given above 
for the path resistance method; however, the argument manipulates the x T Kx term. More specifically, the 
difference across each edge ( i,j ) € K gets written as a telescoping sum of the differences along y t] . The 
Cauchy-Schwarz inequality is applied, and the resulting sum is reorganized in terms of edges in L. 

In the unit priority/unit length case, the two methods can be thought of in terms of the following electrical 
analogies: the path resistance method partitions the conductance (i.e., the reciprocal of the resistance) of 
an edges e € E(G) equally among its incident paths. The resulting paths have path conductances; we find 
a multiplier (r max ) sufficient to increase the conductance of every path so that it is at least as big as the 
conductance of the corresponding clique edge. The edge length method also deals with path conductances, 
but mstead assigns each path a demand equal to the path’s length; this insures that each path has unit 
conductance to support a clique edge. We find a multiplier (p max ) sufficient to make each edge’s conductance 

at least as large as the sum of the demands of its incident paths. In both methods, the multiplier shows up 
in the bound. 

The two methods are duals of each other in the following sense: 

Theorem 4.3. Given a path embedding , the best lower bound for the path resistance method taken over 
all allowed priority assignments and the best lower bound for the edge length method taken over all allowed 
length assignments are the same. 

Proof To show this, we will use a representation presented by Kahale in [Kah96], He defines an 
embedding matrix T as follows: Each row of Y represents a path between a pair of distinct vertices. A row 
entry is 1 if the corresponding edge is in the path and zero otherwise (this is easily generalized to fractional 
paths). It is easy to see that multiplying a vector of edge lengths by T computes the length of each path, and 
that multiplying the result by T r sums the lengths of the incident paths for each edge. Thus multiplying an 
edge vector by r T r and dividing the result for each edge by the original length and taking the maximum is 
a way to compute the value p max used in the edge length lower bound. P T r is clearly a nonnegative matrix, 
and Kahale shows, using the properties of such matrices (see, e.g., [Min88]), that the best lower bound for 
this method is the largest eigenvalue of this matrix. 

The path resistance lower bound can be computed by multiplying a vector of priorities by IT t , then 
dividing the result termwise by the original priorities and setting r max to the maximum resulting value. The 
same arguments about nonnegative matrices apply. However, rr r and r r T have the same set of nonzero 
eigenvalues; this is easily seen by showing that for every eigenvector u of a nonzero eigenvalue of rr T , r T u 
is an eigenvector of r r T with nonzero eigenvalue; likewise, if u is an eigenvector of a nonzero eigenvalue of 


7 



r T r, Tu is an eigenvector of IT T with nonzero eigenvalue. Thus the best lower bounds produced are the 
same. □ 

4.2. Star Embeddings. Embeddings other than clique embeddings are useful. In particular, we embed 
stars into graphs with zero Dirichlet boundary conditions. Such boundary conditions specify a set of vertices 
whose values are zero. Let G be a graph with such a zero boundary, and let L be its Laplacian. Let X be 
the set of non-zero vectors consistent with the boundary restriction. We are interested in that value of the 

following quantity: 

x t Lx 
A = min — = — 
xe* x J x 


To find a lower bound, we embed a star into G such that every boundary vertex is mapped to the center 
of the star, and every other vertex is mapped to a distinct leaf. Thus we specify a path in G from each 
non-boundary vertex to a boundary vertex. Priorities, congestions, resistances, and r max are all defined as 
before. The following theorem applies: 

Theorem 4.4. 


max 

For brevity, we omit the proof, which is almost the same as the proof in the clique case. For readers wishing 
to construct the proof themselves, we note that for the star, the Rayleigh quotient of any nonzero vector 
with zero at the center is exactly 1. 

Matrices for the zero Dirichlet boundary case often represent boundary vertices implicitly. That is, the 
rows and columns of boundary vertices are deleted from the Laplacian of the full graph G. In that case, 
edges between non-boundary and boundary vertices are not explicitly represented by off-diagonal entries in 
the resulting matrix. Instead, they show up as surpluses in the diagonal entries for the non-boundary end 

of the edge. 

4.3. Paths with Masses. To simplify the computation of quantities such as congestions for embed- 
dings into trees in the unit priority case, we introduce the idea of mass. Consider a tree T and a umt priority 
embedding of a graph H into into T (note that this embedding T is unique). Using the definitions from 
Section 3.3 above, we can compute the critical path. We distinguish two types of vertices that do not lie on 
the critical path: those not in any path containing a critical path edge (noncontributing vertices), and those 
that are in some such path (contributing vertices). For each contributing vertex v, there is a first critical 
path vertex v c on any path that includes both v and a critical path edge. Note that since we are working 
with a tree, v c is the same for every such path. We will call u c the connection point for u. Assume that 
the critical path has k vertices indexed from 1 to k in order. We assign masses to the critical path vertices 
as follows: Vertex v { has a mass m t equal to 1 plus the number of contributing vertices that have v t as a 

connection point. 

We refer to the graph P consisting of a path on k vertices plus the masses for the critical path as the 
path-with-masses model. It has an embedding of a subgraph of H consisting of a subset of T that includes 
paths with both endpoints on the critical path. Congestions are computed as before, except that now each 
path contributes an amount equal to the product of the masses of its endpoints to the edge congestions. It 
is easy to see that the path congestion of P in this model is equal to the congestion on the critical path of 

the original embedding. 
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We can also use the path- with- masses model in computing upper bounds. Let y be a vector of length 
k. Let L(P) be the Laplacian of P, and let M be the diagonal matrix with entry [M]u = m t . Let L be 
the Laplacian of T, the tree we are working with. Construct a vector x on n vertices as follows: For the z th 
vertex on the critical path, set the corresponding entry of x to For each noncontributing vertex, set the 
corresponding entry of x to 0. For each contributing vertex, set the corresponding entry of x to the value of 
y for its connection point. Recall Laplacian Property (3.2) from Section 3.2, and note that the only edges 
of T with nonzero differences across them lie on the critical path. Note also that the number of vertices in 
T assigned the value y % is equal to m*. It is therefore easy to see that 

x T Lx __ y T L(P)y 
x T x y r My 

In the star case with Dirichlet boundary, it is easy to see that if we enforce the boundary condition in y, 
then it is enforced in x, which implies that y y T^ y is an upper bound on A as defined above. For clique 
embeddings, note that enforcing the condition y r A/l = 0 is equivalent to enforcing the condition x T l = 0; 
in this case, ~ r|f y y is an upper bound on X 2 . 

5. Lower Bounds for Trees. In this section we show that, for an arbitrary unweighted tree T with 
diameter diam(T), assigning all paths priority 1 in the clique embedding gives a lower bound within a 
logdiam(T) factor of \ 2 . The argument works by first showing the bound holds for a tree with a single zero 
boundary vertex, then applying this fact along the critical path to handle the clique-embedding case. 

We start with a Dirichlet boundary case, where A is defined as in Section 4.2. Let k be the length of the 
critical path. 

Theorem 5.1. For any unweighted tree with a single zero boundary vertex , let r max be the maximum 
path resistance for the star embedding in the unit priority case. Then 

5£M!!,n (A| . 

f'max 


Proof. Note that if the boundary vertex separates the tree into multiple components, the components 
can be considered separately. We need only consider the component with the critical path. 

Let L be the Laplacian of the tree. For any such tree and boundary, we construct a nonzero vector x 
consistent with the boundary condition such that 

(5 x T Lx ^ c max (log k, 1) 

^ x r max 

where c is a constant independent of the choice of tree and boundary. By the definition of A, the theorem 
then holds. The construction of x is done with respect to the path-with-masses model of the critical path 
as described in Section 4.3. 

Let P be the critical path. By assumption it has k edges. We index its vertices as follows: vertex 0 
is the boundary vertex. Vertex i is the vertex at distance i from the boundary. The maximum index is k. 
Mass m. is assigned to v, as specified in Section 4.3 to account for paths that start off the critical path, but 
intersect with it. 

Because we are dealing with the critical path plus a zero boundary, we can express the value as 

follows: Let Si = x t - 1 . Then x* = Ylj=i This allows us to rewrite some of the quantities from 

Inequality (5.1) as follows: 


(5.2) 


1 = 1 = 1 

r ' max E?=i E}=« to; E.*=i * to. ’ 
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and 


x x T Lx ELi % 

^ ^ xTx r \ 2 

Ei=i m i \£j=i s i) 

We specify x by setting <5* = t“*. For fc = 1 it is clear that the quantities on the right-hand sides of (5.2) 
and (5.3) are equal and the theorem holds. For k > 1, we have 

]T<5, 2 = X/ -1 = 0(]ogk). 


Thus if we show that the quantities in the denominators of the right-hand sides of (5.2) and (5.3) are within 
a constant factor, the theorem holds. Note that the two denominators can be thought of as expressions in 
the mi’s, so it will suffice to show that the corresponding coefficients for each m t are within a constant of 
each other. Thus, we want to show that i is within a constant factor of J” 5 ) • For t = 1 these are 

both 1. For i > 1, the summation can be bounded using integration techniques; in this case the sum lies 
between u and 2fi. Thus the square of the sum lies between i and 4t, which gives the desired result. □ 

It is possible to find examples for which this gap is tight. Consider a tree consisting of a path graph 
connecting the center vertices of a series of stars. In particular, let the path have k + 1 vertices numbered 
from 0 to k. Vertex v r (0 < i < k) on the path is the center of a star with \k 2 /i 2 } — 1 leaves. For simplicity 
of reference, we refer to this tree as the bad tree for parameter k. Vertex v 0 is the zero Dirichlet boundary. 

We show that the gap for the bad tree is tight by showing a set of priorities that give a lower bound 
within a constant of the upper bound constructed as per the proof of Theorem 4.4. Let T be a bad tree with 
fe > 1; let L be the Laplacian of T. The path connecting the stars is obviously the critical path. We assume 
that the path vertices are numbered as in the definition of the bad tree. Path edges are numbered by the 

larger of the endpoint indices: i.e., edge i is from Vi~\ to v*. 

Theorem 5.2. Let Vo be the zero boundary vertex for bad tree T with parameter k > 1. Let be the 

unit priority lower bound estimate on A of L given that boundary . Then X — 0 ^ r ^ * 

Proof. We again use the “path with masses” model. It is clear that vertex v t has mass [£ • We note 

that the quantity ^ is always greater than or equal to 1 in the allowed range, so its ceiling is bigger by no 
more than a factor of 2. Because we are interested in bounds, we drop the ceilings; it is easy to verify that 
this does not change the results below by more than a constant factor. 

As in the proof of Theorem 4.4, we can construct a vector x with x 0 = 0 and, for i > 0, Xi = i 2 . Values 
at vertices off the critical path are set as per Section 4.3. Consider the Rayleigh quotient As in the 

previous proof, the numerator is ©(log k). The denominator is 


k 

x T x = ^ m x 

i— 1 



k 

E m . 



2 


As noted in the proof to Theorem 4.4, j~^ is \fi to within a constant factor, so, to within a constant 

factor the denominator is 

k k h 2 k k 2 

= — =^ — = 0(fc 2 logfc) 

i= 1 i=l »=1 

(note that k > 1 by assumption). Thus the upper bound is 0(pr)- 
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For the unit priority lower bound, it is obvious that the critical path runs from the boundary to vertex 
v k . As noted previously, we can write the formula for r max as follows (once again we drop the ceilings without 
changing the result by more than a constant factor): 


k k 

^max ~ ^ ^ i TTli — XI T = e ( fc2 logfc). 

«=i «=i 1 

Thus the upper and lower bounds differ by a log k factor, the value of the Rayleigh quotient numerator. 

We can get a better lower bound and close the gap by assigning priorities to paths. The relative sizes of 
the priorities determine the relative shares of the conductance that the paths get. In the current problem, 
the resistance of the critical path is high because it gets very small shares of the conductances of the edges 
close to the boundary. By increasing the shares of longer paths, we can offset this problem; because the 
number of shorter paths using edge i grows fast as the boundary is approached, the number of paths getting 

an increased share of conductance is relatively small and doesn’t increase the resistance of these shorter 
paths too much. 

We now specify a set of priorities that gives us a lower bound that matches the upper bound. Each path 
that either starts at v, on the critical path, or that has v, as its connection point, gets priority p, = \fi. 
With this change, the congestion of edge i is equal to 


J 3 = 1 

For all i greater than 1, we use integration techniques to get the bounds 


(5.4) 

and 

(5.5) 


k 2 y T -§ ) 

\ y/i VkTl) 

Vv^r s/k) 


For i 1 the upper bound is replaced by 2 k 2 ^2 — J- j ; the lower bound remains unchanged. 

Recall that in the case where path 7 has priority p, we compute the path resistance as follows: 



Plugging in the upper and lower bounds on edge congestion in (5.4) and (5.5) above gives upper and lower 
bounds on the path resistance for a path 7, starting at v 3 . We start with the lower bound; the third line 
below follows from an application of standard integral techniques for bounding sums: 


1 

V 


£cong(e) = — £]Tr§ 

i=l l=i 

>^y(A 

Vj ~{ \yfi \/k+\) 

-l/S)- 


2 / ^\/T + 1 — 2 
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With some straightforward calculations, the reader can verify that the term in parentheses in the last line is 
greater than ^ for all j and k such that k > j. 

The calculation for the upper bound is similar: 



The sum in parentheses is easily shown to be less than 3 for all allowed values of j and k (i.e. , k > j), 

The combination of these bounds implies that, for every path in the embedding that starts at a vertex 
vj, the path resistance is 0(fc 2 ). Note that for any path starting at one of the leaves this value is increased 
by at most 1; thus any path in the embedding has resistance proportional to k 2 . Therefore the lower bound 
on A is ©(^r), which is within a constant of the upper bound. □ 

We now show how to combine the star embedding lower bounds for paths with masses to get a lower 
bound on the uniform priority path resistance lower bound for trees. Consider what happens when we pick 
any point along critical path of the clique embedding and set it to zero. The zero point splits the path in 
two, producing two components with zero Dirichlct boundaries. We can find a lower bound on the smallest 
eigenvalue of each component. 

To minimize subscript conflicts, we introduce the notation r* to stand for r max . We have three r s 
to consider: one from the clique embedding in the original graph, which we denote as r‘ K , and the two 
zero-boundary bounds with respect to our split point. Since we are working with a path, we can think of the 
path laid out with vertices in increasing order from left to right; we therefore denote the two path resistances 
used in the bounds as r* for the left side and r* for the right. 

Recall that the lower bound for the clique embedding is j|jr; we can relate it to the boundary case as 
follows: For tree T with a clique embedding and uniform priorities, let P be the critical path with vertices 
indexed from 1 to k + 1 and let be the path resistance of P. Let l<s<fc+lbe the index of a split 
point. Removing v a separates T into subgraphs; let T\ and T r be the subtrees containing s plus the vertices 
of P with indices less than s and greater than s respectively. Let rf and r* be the maximum path resistances 
for the star embeddings of T\ and T r respectively when v s is a zero boundary point and uniform priorities 
are used (if either of these trees is empty, its maximum path resistance is zero). 

LEMMA 5.3. For the situation described above , 

1 n 

r{+r* < ^ 

Proof. The critical path has k edges and k + 1 vertices, with edge i between vertices i and i + 1. Since 
all of the terms involved are positive, we work with the reciprocals of the quantities in the lemma statement. 
The reciprocal of the left term can be written as follows: 

s — 1 i 

r l*+^r 

t=l j= 1 i=s j=i+ 1 
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We can write the reciprocal of the right expression as 


r * i K i 1 fc+i 

t=l \j = 1 j=i+l 

s— 1 / t 1 


= 53 53 m i’ 

Z=1 Vj=l 


n 


+£ 


f gk 


fc+1 


l 


• m * 
j=*+i 


s— 1 


< 


fc fc-f-1 
*— 1 j~ 1 i=s j=i+l 


To sec the final inequality, recall that n = Rq ■ This proves the lemma. □ 

We can use this to show how to assemble a vector that gives a Rayleigh quotient for the original Laplacian 
that is within a logdiam(T) factor of A 2 . This is shown in the next theorem: 

Theorem 5.4. The uniform priority path resistance method produces a lower bound ff- < \ 2 for the 
Laplacian of any unweighted tree T that is off by a factor that is 0(logdiam(T)). 

Proof. Let L be the Laplacian of T. We consider splitting the critical path at various points. Since rj* 
and r* depend on the split point s, we make that clear by writing them as rf(s) and r r *(.s) respectively in this 
proof. Note that as s increases, r*(s) increases and r*(s) decreases. We also have that r*(l) = r*(fc+l) = 0. 
Let 0(s) — This ratio is unbounded when s — 1 and decreases as s increases, reaching 0 when 

s = k+ 1. We consider two cases: when there exists an s such that 2 > 0(s) > i, and when there is no such 
s. 

In the first case, Lemma 5.3 gives us the following: 


(5.6) 


1 n 

r t(*) + r?(s) < 2 ' 


Theorem 5.1 gives a way to construct the vectors x, and x r so that the respective Rayleigh quotients 
for the left and right sides are O and 0 respe ctively. Recall that these 


vectors are positive. The left Rayleigh quotient has value 


RQi = 

2-/*=i m i x i 

the right has value 


RQ r = 


E fc+1 o 

z=s+i rnixf 


We construct a vector x for the whole tree as follows: assign the values from x, and the negatives of 
the values of x r to the corresponding vertices on the critical path, with = 0. Extend the vector to the 
rest of the tree by assigning vertices off the critical path the values of their connection points as described 
in Section 4.3. Using the Rayleigh quotient for x as an upper bound on A 2 requires x’s entries to sum to 0. 
This is easily accomplished by scaling one of the vectors xj or x r ; this scaling does not affect the respective 
Rayleigh quotients RQ and RQ r . By construction, we have the following (remember that x s = 0): 


t T Lx 


EZl fa - 1 ) 2 + EL fa - * i+ 1 ) 2 

Ei=i rn iX 2 t + Ef= s Vi rn t x'f 
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This is the sum of the numerators of RQx and RQ r divided by the sum of their denominators. It is well 
known that, for positive values a, b , c, and d, 


a + c 
b + d 


< max 



Thus we have that 


(5.7) 


A 2 < max(/?(5i, RQr,)- 


By (5.6), (5.7), and the condition on /3(s), the desired result holds for this case. 

Now assume that there is no s such that 2 > 0(s) > Then there is some s such that 0(s) > 2 and 
0{s + 1) < l By Lemma 5.3, both and r . ( 8 +l) + r . (a+ i) are lower bounds on and hence on 

A 2 . The second of these plus the condition on 0(s + 1) imply that 3r . ( 2 9+T y < %• A similar argument shows 
that 3r ?( a \ is a l so a lower bound. 

We construct a vector x to demonstrate an upper bound. The construction is slightly different this time; 
we construct x, with respect to the boundary being at vertex s + 1 , and x r with respect to the boundary 
being at vertex s. Corresponding values are mapped into x as in the previous case. Once again, we can scale 
one side as necessary to insure that the values in x sum to zero without affecting the left or right Rayleigh 

quotients. 

Note that if it were the case that 


x 7 Lx mim(RQ\) H- num(i7Q r ) 
x T x ~ denom(fiQi) + denom(-RQr) ’ 

the argument as in the previous case would show that ma.x(RQu RQr) would be an upper bound on A 2 . 
This, combined with the lower bounds for this case demonstrated above, would be sufficient to prove that 
the theorem statement holds for this case. However, there is no zero point in this case, and the equality does 
not hold. To solve this problem, we will show that in fact the Rayleigh quotient is no more than a factor of 

two larger than the right hand side of the inequality. 

To see this, note that the denominator of the actual Rayleigh quotient is the same as in the preceding 
expression. The numerator has the following changes: we lose two edges from the ends of xi and x r to 
presumed zero points; this decreases the numerator by x 2 a + x 2 s+1 . We replace these edges by an edge 
between and w. s+ i that contributes (z a - x s+ \ ) 2 < 2x 2 + 2x a+1 . Hence 


x T Lx 


< 


2 (num(RQi) + num(i?Q r )) 
denom(RQi) + denom(RQ r ) 


< 2 max(RQi, RQr), 


and there are upper and lower bounds within a log factor of the diameter. 

This proves that there is a constant such that the theorem holds for the second case, which implies that 
the theorem holds for both cases using the larger of the two constants from the cases. □ 

6. Open Questions. A number of interesting questions remain open. How good in general is the 
bound based on uniform priorities? Is there an easy way to set path priorities that works well in general? 
There is also the question of whether there is an easy way to calculate priorities for trees in both the Dirichlet 
case and the case without boundary conditions that improves the lower bounds. 
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